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ABSTRACT 


A fractal tree of nested viscoelastic boxes is proposed to describe the elastic after-effect in some 
solids that show evidence of a hierarchical structure, such as clay, some rocks, concrete and 
polycrystalline alloys. A functional equation is derived and solved for the Laplace’s transform 
of the creep kernel of the whole hierarchy. Inverting the transform, Hopkinson’s power law is 
obtained as a (fairly) accurate approximation, after an initial time interval. The exponent in the 
power law is obtained as a well-defined function of the branching numbers and scaling 
parameters of the viscoelastic hierarchy. Then we consider a composite solid with both elastic 
(impact) strain response and power law type elastic after-effect, for an applied stress step. The 
stress response to an applied strain is given in terms of Rabotnov’s fractional exponential 
functions. The stretched exponential stress relaxation to an applied strain step is derived as a 
fairly close approximation. According to the present mathematical model, for a given composite 
solid, the stretch parameter of the stretched exponential should be equal to the exponent of the 
power law. The experimental lower bound (0) and upper bound (1) for this parameter are 
derived from the theory. A viscoelastic generalization an equation proposed by Love to describe 
the longitudinal vibrations of bars of uniform cross-section and including the effects of lateral 
inertia, is briefly studied to describe mechanical vibrations in power-law solid bars. 


Key Words: Mechanical vibrations of solids, elastic after-effect, fractal hierarchy, power law 
strain, stretched exponential stress relaxation, functional equations, Rabotnov’s fractional 
exponential function, Love’s equation for viscoelastic beams, non-destructive testing and 
diagnosis, resonance experiments. 


INTRODUCTION 


In the framework of continuum mechanics, a (fairly) detailed study of mechanical 
vibrations can be done, using an analytical approach, if the body is homogeneous, 
isotropic and with linear viscoelasticity (linear hereditary properties), and if a simplified 
enough mathematical model is used. 

When the method of separation of variables can be applied to solve the dynamic 
equations of the mathematical model, displacement fields can be expressed as a series in 
terms of space Eigen functions weighted by the corresponding time functions. With this 
construction both, boundary conditions and initial conditions can be satisfied. 


It is an interesting and general result that the time functions T (t) are the solutions of 


linear integral-differential equations like this one (Rabotnov, 1980): 


! This archive is a modern version of the paper. The original article was written with an IBM 
electric machine and the formulae were handwritten. The presentation was done using a 
blackboard in the old way. Transparent sheets (“filminas”) for text projection and formulas were 
available, but of course there was no power point in those times. 
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The frequency @ appears during the process of solving a suitable boundary value 
problem for the space functions. 
The function q(t) stems from the projection of the density of the distributed external 


—00 


forces applied to the body onto an Eigen function and is zero if these forces are absent. 

The kernel L(t) characterizes the hereditary (delayed viscoelastic) properties of the 

material. 

As an example, let us consider the following equation derived by Love to describe the 

longitudinal vibrations of bars of uniform cross-section and including the effects of 

lateral inertia, to describe mechanical vibrations in power-law solid bars (Love, 1927): 
6’s(t,x) 0°s(t,x) Ge \ 6*s(t,x) 

Pee ax? 010) “ax? OY 


The bar material is assumed to be linear elastic, homogeneous and isotropic, being 
s(t,x) the longitudinal displacement (a function of time ¢ and spatial position x along 


=iE. +pX(t,x) (2) 


the bar axis), o the density, E, the Young modulus, v, the Poisson modulus, and 
xX (t,x) gives the distribution of longitudinal external forces. 


The characteristic radius 7 is given by a double integral taken over the bar cross- 
: 1 
section S,,, being A(S,,) its area: i= WS.) | | ( y+z )dy dz (3) 
cs} S.. 


The origin y =0 z =O is located on the bar axis. 


The viscoelastic generalization of equation (2) can be done as follows. 
Let us consider a linear hereditary relation between the longitudinal scalar stress 


a(t, x) and the corresponding longitudinal scalar strain ¢ (t, x) = = s (t, =) : 


Se [e(ox)- { Ler) (4) 


= 


; : 1 f 

Its inverse is: H(os)=p oth)» | K(-ta(es)a (4 b) 
0 00 

Then, a possible viscoelastic generalization of equation (2) is this one: 


0*s (t,x) a a i L(t-t') 6°s(t',x) aro ENGI) iss (5) 


at? an ace : ax? at 
. E, 

In equation (5), c, is the velocity of extensional waves: Cy =,/— (6) 
p 


Now, let us consider a particular case of the constitutive equations (4 a) and (4 b). 

It is known that certain solids, that show evidence of having a hierarchical structure, 
like clays, some rocks, concrete and certain polycrystalline alloys, often present 
hyperbolic delayed strain responses (known as power law creep) to a step of stress, and 
exponentially stretched stress responses to a step of strain, at least very approximately, 
during long periods of time, after a suitable initial time interval. 
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After its discovery in 19th century’, both hyperbolic and stretched exponential 
responses were often found studying the mechanical properties of many bodies, 
including glasses, polymers, rocks, concrete members, polycrystalline alloys and 
samples of biological tissues. 

When a constant stress of magnitude o, is applied at the origin of the time variable, 


being zero before, the strain response é(t) to an applied step of stress has a linear elastic 
component ¢, anda delayed component g(t): é(t) = & +é,(t) 


is ; ; : Mia ts 
IfJ, = a is the elastic compliance modulus, the elastic response is given by €) =J)-o> 
0 


Dp 
The delayed strain seems to follow a hyperbolic law ¢, (t) ~C- a ‘Oo 


Here C is a positive parameter. This power law fits very accurately experimental data 
for times ¢ greater than the characteristic timer, . The parameter p is always between 0 


and 1. 
So, the whole strain response to a step in the stress is given by: 


p 
E(t)=] J, -e{4] Oy 
For a variable applied stress o(t) the delayed strain response of a linear viscoelastic 
power law material is given by: E,(t)= fol —t')-o(t')-dt! 
In this case Q(t) =k (t) is the hereditary kernel in the case of this kind of material. 


If o(t)=0, for t>0 we obtain ¢,(t)= 0, a O(t')-dt' 


Dp 
Then froeyar=cf 2 so the viscoelastic kernel is (with a@=1l-—p and 


o(7,)= 2b") o(r)=a(7.) (=) a 


If o(t)=0, for 0<t<T, and zero otherwise, then the delayed strain response for t > T, 


« nrarffor-rrarcave(z)”(1-(-B)’ 


But in this case, if t — +co the delayed strain response €, (t) > 0 


As, of course, the instantaneous elastic response disappears also, the material behaves 
as a solid. This behaviour justifies the usual denomination of “power law solid” applied 
to this kind of materials. 

Now, notice that the power law response can be writtene,(t)=VJ/(t)-0, 


with J(¢)=J(1)-t?.The compliance J(t)verifies the following functional equation: 
J(ty-x)=M(t,)-F(x) Here M(t,)=1? and F(x)=VJ(1)-x?. 


2 Apparently, the power law creep was discovered by W. Weber in 1836, studying the delayed strain 
response of the silk threads of galvanometers under a constant tensile load. 
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With ¢, fixed, x [0,1] andt=t¢,-x e[0,t,], this functional equation characterizes a 
scale invariant response J(t) in a time interval[0,t,]: M(¢,) gives the amplitude and 


F (x) gives the shape (always the same) of the response. Conversely, the only regular 


solutions of the functional equation of scale invariance are the power functions (Pipkin, 
1972). 
The stretched exponential stress relaxation is given by: 


9, (t)~ o(0)-exp {4} (8) 


For times ¢ greater than7., the stretched exponential fits to experimental data fairly 


accurately. The parameter p also is always between 0 and 1°. 


There is another fact that is closely related with the hyperbolic response. It is the 
relative constancy of the quality factor corresponding to elastic wave dissipation in the 
Earth, for most earth materials and for wave periods comprised between 10° and 10° 
seconds: it lies between 100 and 1000 (Pilant, 1979). 

As will be seen here, a possible interpretation of the hyperbolic or stretched exponential 
behaviour of power-law solids can be obtained assuming that strain and stress responses 
result from the addition of partial strains or stresses that correspond to the degrees of 
freedom in an infinite hierarchical structure. 


The purposes of the present paper are twofold: 


(1)-First, to apply the elementary mechanical theory of composite materials to a suitably 
defined hierarchical structure of elastic elements of different sizes and shapes immersed 
in a viscous matrix. A given degree of freedom will be able to proceed to a relaxation 
process, if the precedent degree of freedom has relaxed before. As consequence, the 
time scale of evolution of a given level is subordinated to the time scale of the precedent 
level: each level responds slower than the precedent level but faster than the next one. 
Under the assumption that the dominant viscous dissipation occurs at the interfaces 
between the elastic elements and the viscous matrix, we obtain the power law creep that 
describes the delayed strain response of the above-mentioned composites. 

In order to do that, we derive a functional equation in Laplace domain, whose solution 
in time domain corresponds to power law creep. This derivation follows an approach 
closely related with Liu’s explanation of constant phase angle impedances and 1/f 
noises in rough electrode interfaces (Liu, 1985). 

Then we show that the stretched exponential can be considered as a fairly good 
approximation for stress relaxation in a composite with delayed power law and non-zero 
impact response to a stress step, after an initial time interval and during a significant 
part (but not towards the end) of stress relaxation. We derive the same value for the 
parameter of the hyperbolic response and for the parameter of the stretched exponential, 
as well-defined function of the basic parameters of the viscoelastic hierarchy. 


3 The stretched exponential seems to have been discovered by R. Kohlrausch in 1847 while studying 
the relaxation of shear stress under conditions of constant shear strain, also in the silk threads of 
galvanometers. 
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(2)-Second, to apply equation (5) to resonance experiments done with bars of power- 
law solids. 

We separate variables in the case of a bar with free ends, assumed supported by its 
midpoint. 

We derive approximate formulae for resonance frequencies. 

Mechanical resonance or free vibrations produced in bodies of suitable shapes, sizes and 
material allows very accurate measurements of Young’s elastic moduli. This is the case 
of concrete, glass, ceramics and metallic alloys. Longitudinal, flexural and torsion 
vibrations may be excited and measured in order to estimate the elasticity moduli (see, 
for example, Trianon and Griego, 1983 and the bibliography suggested in this report). 
The following formula is used to estimate Young’s modulus from the measurement of 
the fundamental frequency of resonance f, of longitudinal vibrations, the density p of 


the material and a correction factor C, : b= pC, i (9) 


For power-law solids, we suggest a way to incorporate the effect of the internal friction 
in the above mentioned correction factor. 


MATHEMATICAL MODELING 


(A)-Power-Law Solids 


Let us consider a biphasic composite material. One phase is a linear viscous matrix. The 
other phase is a collection of linear elastic aggregates composed by many elements of 
different sizes and shapes. 

We apply a normal stress o, =oon the opposite faces of a volume element. Let 


€, = €be the resulting strain in the direction of the applied stress. (Figure 1) 


Figure 1.Element of volume of a biphasic composite with hierarchical structure. 


Working in the framework of the elementary theory of composite materials, we suppose 
that the connection between stress and strain can be obtained in one dimension, for this 
composite, using a sequence of nested viscoelastic boxes (Figure 2). 
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Figure 2.Viscoelastic behavior of the composite represented by a set of nested viscoelastic boxes. 


h-+u 


The viscoelastic box of order n is composed of Nnt+lelements in parallel: one linear 
viscous element with damping parameter 7,and Nn complex viscoelastic elements. 


Each of these complex elements is the in-series combination of a linear elastic element 
with compliance J, and a viscoelastic box of order n+1. 


Then, ifo, (t)is the stress applied to the viscoelastic box of order n at time, and if ¢, (t) 
is its corresponding strain response, we have: 


t)=[O,(e-2)-0,(¢)-ae (10) 
We assume that the stress and strain are always zero when t<0. The creep kernel of the 
viscoelastic boxes of the level n in the hierarchical viscoelastic structure is Q,(¢). 


The viscoelastic box of order 1, which contains all the other boxes nested in its interior, 
is supposed to be in series with a single elastic element of compliance J,. This element 


gives the impact or immediate elastic response of the considered element of volume of 
the composite material. If ¢', (t) is the strain produced in the elastic element of 


compliance J,,, then: é,(t)=e', (t)+é,,,(t) (11) 


n? 


Bute’, (1)=J,-0,,(¢), sothat: — 4,()=J,-0,,0)+[ O,s-2)-o,(¢)-de’ (12) 
But in parallel to the combination of elastic compliance with viscoelastic box of order n 


there is a viscous element. 


The stress acting through it verifies: o,,(t)=7, : Ke, (t) (13) 


Ifg,,is the volume fraction of this viscous element, referred to the whole volume of box 


number n, and if g,is the volume fraction of the in series combination of an elastic 


element with box number n+1, from the elementary theory of composite materials we 
have, applying the popular rule of mixtures: 


On()= Ny Py Fnal)+ Pry" Fall) (14) 
The volume fractions verify, for every n: NP, +Qry =1 (15) 


The possibility of applying the rule of mixtures, at least as a convenient approximation, 
can be justified as follows. If A, is a cross sectional area corresponding to the boxes of 
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level n in the viscoelastic hierarchy and A,,is the area corresponding to the dashpot, 


A Ay, 
n a . Ons (v)+ he ° o,,,(t) 


then: 4,-o,(t)=N,-4,,,-0,,,(t)+4,,°¢,,(¢) so o,(t)=N, - 


n+l 
n n 


A nv 
As4,=N,-A,,,+4,,,80 we have N,-—4+4+—=1 

A, A, 
However, the average fraction of the area occupied by a certain component in a set of 
parallel cross sections of a composite body can be used as an estimation of the volume 
fraction of the component in the body (Delesse’s rule; Weibel and Elias, 1967). 


As consequence of this connection between area averages and volumes of a given phase 


n+l 


in a body of composite material, if is interpreted as an estimation of the average 


n 


fraction of parallel cross sectional areas of viscoelastic boxes of level n+1 inside a box 


ntl 


of level, taking into account Delesse’s rule, we have @, = 


nN 
nv 


A 


n 


Of course both stresses and strains must be 


An analogous derivation gives 9, = 
interpreted also as suitable averages. 


(A.1)-Derivation of a functional equation for the creep kernel Q, (t) 


Taking Laplace’s transforms, and defining f (s)= [ ea ee (t): dt, it follows from (11), 


(12), (13) and (14): é (s)=0,(s)-6,(s) (16) 
G,,(s)=s-7, -&,(s) (17) 
E,(s)=J, -Fuls)+ Quuls)-Fa(s) (18) 
G,(s)=N,-Q,-Fruls)+ Py“ Fl) (19) 
Eliminating ,(s), &,,,(s)andG,,(s) between (7), (8), (9) and (10) we obtain: 
= 1 
O,(s)= oe. (20) 
OS FeO) 


Now, if we assume an infinite hierarchy of nested viscoelastic boxes, then from (20) we 
obtain the following continued fraction: 


=~ 1 
Q(s)= a (21) 
; vA 
Py Mi S+ 1 
ie 
Q,, My -S+ N° 
2v 2 oi + 
Now, let us suppose that: VN, =N J,=J/ 90,=9 @,,=@,, for every n. 
Then (15) reduces to N-9+9,=1 (22) 


If the viscous dissipation of elastic energy is produced mainly in the region of matrix 
surrounding each aggregate-matrix interface, we can expect that as the aggregate 
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becomes smaller, the viscosity becomes larger (the other conditions remaining 
unchanged). 
So, introducing a parameter ( between 0 and 1, we assume that the viscosity 


coefficients of the different levels of the viscoelastic hierarchy verify: 
n-1 
n,-(¢°) =n (23) 
1 
So, as the volume fraction @ is less than 1, 7, grows scaled by —; from each level to 
2 


the following level of the hierarchy. Then (21) reduces to: 
1 


Q,(s)= Ne (24) 
PMS + ] 
os 
—, th “St+ ae 
a J+ 7 a 
2, 1S S+ 
(0°) J+ 


Taking (14) into account, and splitting the volume fraction g = vy’ -y'”,, it is possible to 
re-cast the continued fraction (24) as follows: 


= 1 

Q,(s)= N-o? 
PM St J 1 
at 


(25) 


N-o'? 
PM ST J 1 


+ 
(0° J QO, St... 
This expression is convergent because N-g@ <1 (Wall, 1972). 


Pp 


For a given set of parameters of the viscoelastic hierarchy, (16) defines a function of s 
for positive values of the Laplace’s variable. 

If we substitute in (25) s by g’-s we obtain, after some operations, the following 
functional equation: 


(op. 
o a s) _ ! (26) 
gr? 
gp? . T ‘Ss + ~ Fa) 
1+ O's) 
J 
The characteristic time 7 is defined by the equation T=9,:7,-J (27) 


Equation (26) can be solved numerically. However, here we want to proceed analytically 
as far as possible, so we are going to derive an asymptotic functional equation whose 
analytic solution is known. As also shown in the mathematical Appendix, assuming that 
n-g' >1, when sJ0 then O,(s)T oo and s-O,(s)/0 As consequence, when 7'-s is 
small enough relative to 1 we will have £als) small relative to 1 and as) large 
relative to 1, so (26) may be approximated by: 
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a = 
Ole" -s)= Gils) (28) 
Regular solutions of this functional equation are of the following type: 
= A 
Q,(s)= P (29) 
The parameter p is given as the following function of some parameters of the viscoelastic 
I = 
hierarchy: Sy EL (30) 


<7 
In} — 
QP 


As will be shown later, p is a parameter comprised always between 0 and 1. 


Inverting the Laplace transform, we obtain the desired delayed strain kernel that gives a 
power law strain response to a step of stress input: 


t a 
O(t)= ant) (31) 
Here a@=1-p (also comprised between 0 and 1), and ¢, is an arbitrary but positive 
reference time. If we taket, =7.., we can expect that (31) approximates the viscoelastic 
strain kernel fort > 7. However, C remains undetermined, as well as Q, (T,) which is 


proportional toC . 


(A.2)-Derivation of the stretched exponential stress relaxation 
If a linear viscoelastic material has both instantaneous (elastic) ¢,(t) and delayed 


é,(t)strain response to an applied stress o(t)that is zero for negative times, the total 
strain ¢,(t)+ ¢,(¢) isgivenby: — e(t) =J, (o(e)+ fix (#1) o(t') dt ) (32) 


If the case of the composite material studied in this work, the creep kernel K (t) may be 
identified with Q,(t), the kernel of the box in level one in the viscoelastic hierarchy. 
Inverting (32) and introducing an elastic impact modulus £, such that J,-£, =1 and a 


stress relaxation kernel L(t): o(t)=E£, (2()- fiz(e-#/) é(t') ar’ (33) 


Taking Laplace’s transforms in (32) and (33), it is possible to relate the transform L (s) 
of the stress relaxation kernel with the transform K(s) of the creep kernel: 


ye) 


“T+ K(s) so 


: = A 
Now we substitute the Laplace’s transform Q,(s)=— of the power law kernel (see 
s 


7) 
equation (29)) in place of K (s) in (34), wit = [z) , and we obtain: 
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£(s) Rica (35) 


1 
Pp 7 i, a 
From 2 = [z) and from (34) it follows that: y= (=) pe { Tp) i (36) 


0 é 


Inverting the transform (35): 


ott) Eacrale) 


As usual, T(z) represent the Gamma function. An equivalent expression was called 


fractional exponential function by Rabotnov, already in 1948, since for p =1 reduces to 


t 
agit 
the decaying ordinary exponential. In our case, when p=1 we obtain: [z) -e CZ 


é 


The function yO uy "is approximated by (Rabotnov, 1980): 


L()=2| fr, i .p-(2- py cx] -wa-ay (22) | (38) 


an : 
The relaxation of stress at a constant strain ¢, for a viscoelastic body is given by 
o(t)=M(t)-¢, with, by definition: 


_ = f '\, ' 
M (t)=E, (1 Jz) at’ (39) 
Taking into account (41), from (40) it follows that the stress relaxes according to a 
stretched exponential: o(t) =M (t) "Ey =O, elt) (40) 
In this equationo, = M, - 


From equations (38) and (40) we obtain the characteristic stress relaxation time 7, 


given as function of the characteristic time 7, of the hyperbolic strain response: 


a ae (41) 
em TCD 


(B)-Mechanical vibrations and resonance from the generalized Love’s 
equation for power-law solids 
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Let us consider the viscoelastic generalization of Love’s equation in absence of 
distributed external forces, so X (t,x) is always zero. 


Then: 
O's (t,x O’s(t,x) + 0’s(t', O's(t,x 
tng funn ented — 


Now, if the material behaves as a power-law solid, the kernel L(t) is given by a 


fractional exponential function (37) that can be approximated by equation (38). 

Let us consider a bar of length¢, and put the origin of coordinates at one end, 
soO<sx<f. 

If the bar has free ends, the longitudinal stresses verify o(t,0) = o(t,0) =0 for each 


instant of time. 


From equation (4 b) it follows that the longitudinal strains verify Neumann boundary 
O 


conditions: é(t,0) = 5, 5\9) =0 é(t,£) 7 Xs(1,0) =0 for every instant of time. 
x x 
Separating variables by introducing the ansatz s (os) =T (¢) S (x) in equation (42), from 


the boundary conditions we obtain the normalized Eigen functions: 


S(x)=S, (x)= [2 cos(k,, x) k,, = ae m =1,2,3,... (43) 


To each Eigen function it correspond a time function T (t)=T, (t) that verifies the 


m 


integral-differential equation: 


2 t 
a (a8, |Tal)= | ern, (¢)a"] =0 (44 a) 
In this equation: Oom = — (44 b) 


If T,,(t) is zero ift <0: 
[ £(-e')7, (e)ar=[o(- t') (e)ar=f2) (t—1")dt' 


Then the equation (44 a) can be ies by the Laplace’s ‘cneioan in terms of the initial 
d 
conditions 7, (0° and —7 (0°) for7, (t). 
m ( ) dt m ( ) m ( ) 


If the internal friction of a power-law solid is low enough, the following approximate 
asymptotic formula for the fundamental frequency of free oscillation can be derived 


from (44 a), witha(p)=p-(2-p)”: 


é 


l-p P 
+00 T. t 
O =, 1] S [| a(p) oo etn 22 foto t') dt' (45) 
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DISCUSION 


(D.1)- The parameter p of the power law strain response 


and a are 
Q 


From (32) we conclude that p must be positive, because both & = 7 
"Q 

greater than 1. 

But from both, experimental results and thermodynamic requirements related with 

elastic energy dissipation, we know that p must be less than 1. 


This is equivalent to the restrictionv > 1. Taking into account the definition of v given 
in equation (29), as well as the inequality V-~ <1, we see thatv >1 is verified if and 


only ifN-g'? >1. 
The physical origin of this last inequality and the physical origin of the scale law for the 
dashpot viscous moduli 7, (g? i =7,, can be found examining the geometric scaling 


laws of the viscoelastic hierarchy. 

Let us consider an element of volume V,, in the composite material corresponding to a 
combination in series of an elastic element and a whole viscoelastic box of level n. 

Let us introduce a characteristic length/, corresponding to volume V, and a 


n 


characteristic area S, ~/° of the interfaces between the volume element and its viscous 


. . 1 2 
matrix environment. We assume that/,, ~ yp” ‘1,80 S...® yp” -S,, and V, 


n+l = Q : E, 
Furthermore, we suppose that the force F, taken by the viscous matrix in the volume 
0é, 


Ot 
As consequence, the force per unit volume o,,, is given by: 


n F S, O€, -; es S, 0€, 
On, = = i . . = . ese era ats = 
V, V, «Ot V, ot 


We suppose now that the viscous stress ois related with the density of volumetric 


element is given by F, = y-S,,- 


S O& 
viscous forces by a constant factor: of,,=y-:o,, So thato,,=7,- ap , with 
t 


I n-l S 
n=(0" | “mand 7, = 7: ee 
V, 


As consequence, we obtain the scale law of viscous dissipation 77, (po? as =7, with 
1 ; 
B= 3 According to this model, the parameter of the power law and of the stretched 


ner) 


exponential verifies p = ———< 


nf e | 
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, 2 
NS ia a N-0 =(N-9'?) is greater 


than 1: this happens if the total area related with viscous dissipation increases from one 
level of the hierarchy to the next one. 


Now v is positive and p is less than one if 


(D.2)-The relaxation times in a self-similar viscoelastic hierarchy 

The self-similar nature in physical space of the viscoelastic hierarchy proposed here 
gives rise to power law responses, and power law responses point to scale invariance in 
time. Scale invariance, at its turn, implies a manifold of time scales in the delayed strain 


) eal) 
= J(1)-t? - h that the local t Ea 
responses €, (t) ( ) t" -Oy, SUC a € local time scale lde,,(t)/ail 
le, (0) t 


lde,(t)/di] p 
The characteristic times 7, do not appear in the local time scales. 


of the response 


at the instant tis always proportional tor: 


In the case of stretched exponential stress relaxation g (t)~ a(0}e| -[ Vr 


o..\t pe 
time scale is: mi -(%) t 
|do, (t)/dt| t 
It grows liket'’, slower than the growth of the local time scale of the delayed strain 
response. But now the characteristic time 7, appears in the formula for the local time 


P 
} |i local 


scale: the stretched exponential does not have scale invariance in time. 
(D.3)-Correction factor to estimate Young’s modulus 


When the masses of transducers used to excite and register longitudinal vibrations in a 
bar are negligible, and internal friction is negligible, formula (45) reduces to 


2 2 1 or k;, 
hi Joi Cotrenra (46) 


m 


From equation (6), it follows that Young’s modulus Young’s modulus £, is related to 
the density and the velocity of extensional waves: ES pe. 


If Poisson’s modulus is known for the material of the bar, Young’s modulus can be 
estimated from the fundamental frequency using equations (9) and (46): 


2: 
E,=p C, f2 c,-4e[ rele) | 


If internal friction must be considered, from formula (45) we derive E, =p C, f° , but 


2 2 
now with C, = 7 ze [1 (% i) | 
(1-f, L(t')cos(2fy, t') ar’ é 


Through f,,, Young’s modulus appears in the integrand. 
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So the formula above could be applied to relate Young’s modulus £, with the measured 


fundamental frequency in resonance only if a good enough preliminary estimation of 
the modulus £, is available to calculate an approximation to f,. 


It would seem that the analytical approach, as applied in this work, should be 
abandoned. 


CONCLUSIONS 


(a) An analytical formula for the power law parameter as function of the branching 
number and volume fractions of the hierarchical composite material was derived 
for a simplified example. The stretched exponential was obtained as an 
approximation to stress relaxation when strain creeps according to a power-law. 
In this approximation the parameter of the stretched exponential turned to be the 
same as the power law parameter. 

(b) The functional equation method does not allow us to derive an analytical 
formula to relate the constant C that appears in the creep formula 


Pp 

t : : : : : 

&; (t)*C {<| -o, with the parameters of the viscoelastic hierarchy. This 
poses a new research problem. 

(c) A formula was obtained that relates the characteristic times 7, and 7, with the 


mechanical parameters of the composite material. 
(d) Due to its scale invariance, the hyperbolic strain response has an infinite and 


, t walt 
unbounded set of local time scales—. The characteristic time 7, must not be 
P 


considered as a time scale of the hyperbolic response. It gives a measure of the 
numerical order of the interval of time (since the application of a step in stress) 
after which the power law can be applied. 

(e) The characteristic time 7, gives a measure of the numerical order of the interval 


of time (since the application of a step in strain) after which the stretched 
exponential law for stress relaxation can be applied. The stretched exponential 


Dp 
also has an infinite and unbounded set of local time scales of -t Contrary 


to the hyperbolic response case, now 7. appears in formula for the local time 


scale of stress relaxation. 

(f) In principle the equation (12) for the Laplace’s transform of the creep kernel 
allows to take into account random branching numbers as well as random elastic 
and viscous moduli. The resulting equation could be approximately solved 
applying a method developed by Kaplan and Gray to study the effect of disorder 
on a fractal model for the alternating current response of rough interfaces 
(Kaplan and Gray, 1985). 

(g) Finite viscoelastic hierarchies show, besides the high frequency cut off that 
corresponds to the beginning of the deformation process, a low frequency cut off 
for long times of deformation. In the model based in a set of nested boxes, the 
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low frequency cut off behaviour depends of the termination (elastic elements or 
dashpots) of the last box in the hierarchy. 

(h) From the difficulties found with the analytical approach to mechanical vibrations 
intended here, it seems that a fully numerical approach to the generalization of 
Love’s equation for power-law solids must be developed. 
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Dr. S. H. Liu February 26, 1988 
Solid State Division 


Oak Ridge National Laboratory 
Oak Ridge 
Tennessee 37831 USA 


Dear Dr. Liu!: 


Thank you for your kind letter of February 4, 1988; for your letter of November 13, 
1986, and for the reprints on your theory about the ac response of rough interfaces! 

It seems to me that my letter of November 24, 1987 is lost. So I will reproduce here 
what I wrote there. 


Certain of your results about rough interfaces were applied here, in Uruguay, to improve 
the design of the tip of a stimulating electrode of complex shape which is now being 
used for medical purposes. 

While I was studying the reprints I perceived a close connection between your 
explanation of constant phase angle impedance and a possible explanation of 
Hopkinson’s power law for viscoelastic solids. 


Consider a biphasic composite material. 

One phase is (mainly) a linear viscous matrix and the other phase is (mainly) a linear 
elastic aggregate. The aggregate is composed of many elements of different sizes and 
(perhaps) shapes. 


Let us suppose that we apply a stress o, =o on opposite faces of a volume element of 


the material. 


! The original letter was hand-written. The pictures that appear here are computer scans of the 
original drawings. The fractal model explained in this letter was previously part of the presentation 
by R. Suarez Antola, “Vibrations of viscoelastic, power-law solids”, in the First Regional Seminar 
of Latinoamerican and Caribbean on Vibration Analysis and Diagnostics, Bariloche’s Atomic 
Center, Bariloche, Argentina, September 14 to 18, 1987. 


Let ¢, =€ be the resulting strain in the direction of applied stress. 


Working in the frame of elementary theory of composite materials, we suppose that the 
connection between stress and strain in one dimension can be obtained using a sequence 
of nested viscoelastic “boxes” as shown in the following picture: 


hs 
Rhye | 
] foe fe] wl 
| 
te, 


Viscoelastic “box” of order 1 is composed by N, +1 elements in parallel: one linear 


ey 


viscous element with “viscosity” 7, and N,, complex viscoelastic elements. 


Each of these elements is a combination, in series, of a linear elastic element with 
elastic compliance J, and a viscoelastic “box” of order n+1. 


If o,,(t) is stress applied to viscoelastic box of order n at time ¢ and if ¢,(t)is its 
corresponding strain, we have: é,(t) = ie Q,(t-t)o,(t)dt [1] 
If g,,is the volume fraction (referred to the whole volume of box number 7) of the 


viscous element in parallel, and if @, is the volume fraction of the in-series combination 


of an elastic element and a box of order n+1, we suppose that: 


0, (t)=N, P,, Ons (t)+9,, ony (1) [2] 
Here o,,, (t) is stress applied to the viscous element, and  =N,9,+9,,=1 [3] 
Besides: Catal ican ()+[ 0... (t—r)o,,,(7)dt [4] 
And ,, (t)=n, Ke, (1) [5] 


Taking Laplace’s transform in equations [1], [2], [4] and [5], after some work we obtain 
(sis Laplace’s variable): 


é,(s)=9,(s) &,(s) ; O,(s)= NG. [6] 
Pry er FAO.) 


In order to simplify the analysis, from now on we shall take V,=N, J, =J , 9,=9, 


?,, =@, for every n. 


But if viscous dissipation of elastic energy is related mainly with the region of matrix 
surrounding each aggregate-matrix interface, we expect that, as aggregate becomes 
smaller, viscosity becomes larger (the other things remaining unchanged). 


If we suppose that: n= [+] n (22 152,355.2) [7] 
Q 


It is possible to show, beginning with equation [6] and taking n=1: 


= 1 1 


O,(s)= [8] 
1+ O,(s) 
i 
Now, let 7, = 9,7 J be a characteristic time of the material. 
Then, it 7.s<1, az) <1 and ae) >>1, we obtain the same basic functional 
equation of your theory of a.c. response of rough interfaces: 
- fo 
s)=—QO(s 9 
O.(9 ) ono) [9] 
; ; = A 
Proceeding as you have done, we find the solution O,(s)=— [10] 
S 
Where p=li- log. N 


ve,(2) [11] 
log,| — 
QP 


As, by equation [3], 0< Ng@<1, we always have N< = . Then O0<p<l Besides, if 
Q 


N>1 then p<l. 


0 if t<0 


we have: 
0, if t> +H 


Applying a Heavside’s step stress history a(i)=| 


= = A 
&(s)=0,(s) 7a = fay % 


Inverting the transform we find the well-known version of Hopkinson’s power law: 
A : : 
é(t) =| ——— |t” 0, _ where T(z) is Euler’s gamma function, and 0< p<1. 
Tr ( pt 1) 


Experimental values of p are in fact bounded between 0 and 1. But now p, instead of 


being an empirical parameter is expressed, by equation [11], as a definite function of 


branching number WN and scaling parameter a 2 


1-9, 
It seems to me that this approach can be generalized in at least two complementary 
directions: 


First, it can be randomized as your theory was in the paper of Kaplan and Gray. 


Second, it can be extended to a nested hierarchy of viscoelastic fractals with varying 
fractal parameters, so that each one has its own p. In this case, as time grows (or as 


frequency lowers), “equivalent p” will be modified as is actually seen during certain 


experiments with polymers and (perhaps) with concrete? 

As far as I know, power-law models are used in experimental work with viscoelastic 
materials because they are simple and fairly accurate, and in theoretical work because 
they allow us to use the well-developed mathematical framework of Abel’s and 
fractional exponential operators. 

But if this approach is correct, it seems that it can be used to design viscoelastic 
materials with pre-stablished p values and to explain experimental variations in the 


exponent of Hopkinson’s power-law in terms of structural variations in the material. 
Perhaps this kind of approach can be applied to polycrystalline materials under suitable 
conditions, assigning elastic properties to the bulk of the grains and viscous properties 
to the interfaces between grains. 

I must confess that in my personal situation I have certain difficulties to get information 
about what is happening in outside world. So, perhaps someone else has done some 
work in this same direction and I don’t know. 

Anyhow, I can’t develop the experimental aspects of this approach in my country. 

If it is possible, I would like to know if you have published further work on your theory 
of rough interfaces and related topics. 

It seems to me that this approach in temporal domain to viscoelastic power-law suggests 


a natural generalization: Substitute stress o for any suitable spatially intensive variable 
(‘“pervariable” in the terminology of systems engineering), substitute strain ¢ for the 
corresponding conjugated (in relation with energy) spatially extensive variable 


; , 0 
(“‘transvariable’’) and substitute the elastic and viscous laws, €=J o and o=7 a é by 


homologous laws in terms of the new variables. In this case a hierarchy of nested boxes 
constructed in the same fashion used to obtain Hopkinson’s power-law, with suitable 
branching numbers and scaling parameters, will end in a relation between transvariable 
and pervariable which is hyperbolical for times long enough or frequencies low enough. 
Perhaps we have here a possible approach to “Joseph effect” (as defined by 
Mandelbrot)? 


Yours truly 


Roberto Suarez-Antola 


’ 


ee 
POST OFFICE BOE = 


OAK RIDGE NATIONAL LABORATORY 
OAK FIDOE. TENNESSEE 3787" 


OPERATED 8y Marre) MATUETTA ENERGY SYSTEMS, INC 
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Dr. Roberto Suarez Antola 

Analisis Fisicomatematico de Sistemas 
Agracitada 334! 

Montevideo, Uruguay 


Dear Dre Suarecs Antola' 


Thank you for your letter of February 26 which allowed me to underatand 
your interesting idea on the viococlantic properties of composite solids. 
I have done a limited amount of Literature search, and as far as I can tell, 
no one has attempted what you have done. Therefore, I strongly encourage 
you to write it up and subett it for publication in the Phyetcal Review. 
Aleo, if you can find the means to support you, 1 think you should present 
your work at the next Materials Research Society meeting to be held in 
Boston later this year. 1 Will send you the information when I receive it. 


As to the work iteelf, I think it makes eminent sense. Just like my own 
paper, it relies on the initial assumption that there is ao hierarchy of struc 
ture. Purists may object to thie assumption, ae they did when my paper was 
first submitted for publication. Now the fractal nature of surfaces has been 
firmly established. There exists some evidence of hierarchical structure in 
composite materials, notably clay and rock. You may cite these as motivation 
for proposing such a mdel. The last paragraph of your letter touches upon 
a problem I have been studying, i.e. the dielectric response of fast ton 
conductores, which conaist of # random Liquid of fone moving in a regular 
crystalline lattice. Again, a hierarchical atructure of the fontc distribu- 
tion will establish the observed power-law frequency dependence, but strong 
evidence for stich a distribution is lacking. 


1 am also enclosing # recent review of our work on the interfacial 
problem. This represents our latest thinking on this subject. 


Sincerely yours, 


SW Kx 


S. H. Liu 
Solid State Division 


SHL:vh 
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Power Law and Stretched Exponential Responses in Composite Solids 
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Keywords: Solid composites; Elastic after-effect; Fractal hierarchy; Power law strain; Stretched 
exponential stress relaxation; Fractance; Mathematical modelling. 


Abstract. Clay, rocks, concrete and other composite solids show evidence of a hierarchical structure. 
A fractal tree of nested viscoelastic boxes is proposed to describe the elastic after-effects in these 
composite solids. A generalized fractal transmission line approach is developed to relate the strain 
and stress responses. Power law for strain, under an applied stress step, is derived. The exponent in 
the power law is obtained as a well-defined function of the branching numbers and scaling parameters 
of the viscoelastic hierarchy. Then, a composite solid with both instantaneous (linear) elastic strain 
response and power law type (linear) elastic after-effect for an applied stress step, is considered. The 
stretched exponential stress relaxation to an applied strain step is derived as an approximation. For the 
same composite solid, the stretch parameter of the stretched exponential and the exponent of the 
power law result to be equal to each other. 


Introduction 


Certain composite solids, that show evidence of having a hierarchical structure, like clays, some 
rocks, concrete and biological materials, in the linear regime, often present hyperbolic delayed strain 
responses (known as power law creep) to a step of stress, and exponentially stretched stress responses 
to a step of strain, during long periods of time, after a suitable initial time interval ({1] to [7]). 

When a constant stress of magnitude o, is applied at the origin of the time variablet, the 


hyperbolic law for delayed strain is ¢,(+)~ «{4} -o, being K a positive parameter. This power law 


é 


fits very accurately experimental data for times ¢ greater than the characteristic time c,. The 
parameter p is always between 0 and 1. 
For a variable applied stress o(r) the delayed strain response of a linear material with power-law 


viscoelasticity, is given by <,(1)=[,O(t-r)-o(")-d'. The viscoelastic kernel is: at)=ote.){ 


a@=l-p Qcz,)= (l-a)-K The stretched exponential stress relaxation after a step of strain is given 
T 


é 


fom 


q 
by caatdveo | For times ¢ greater thanz,, the stretched exponential often fits to 


experimental data fairly accurately. The parameter g also is always between 0 and 1. 


A hierarchical and fractal viscoelastic model, to describe the mechanical responses of certain 
composite solids, was recently proposed [8]. A functional equation for the Laplace’s transform of the 
viscoelastic strain kernel was derived, giving the desired asymptotic power law strain. For solids with 
elastic impact responses a delayed stretched exponential stress relaxations was obtained. However an 
indeterminate constant appears in the solution of the functional equation for the power law strain and 
reappears in the stretched exponential stress relaxation. 

The purposes of the present paper are threefold: 1-To review briefly the above mentioned fractal 
viscoelastic model of composite solid. 2-To apply a generalized fractal transmission line approach 
(initially presented in [8]) instead of the functional equation approach detailed in [8]. 3-To derive the 
power law and stretched exponential responses without the appearance of undetermined parameters. 
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Mathematical Modeling of the Viscoelastic Hierarchy 

Let us consider a biphasic composite material. One phase is a linear viscous matrix. The other phase 

is a collection of linear elastic aggregates composed by many elements of different sizes and shapes. 
We apply a normal stress o, =o on the opposite faces of a volume element. Let ¢, = ebe the 

resulting strain in the direction of the applied stress (Fig.1 left). 


Fig.1: Left: Element of volume of a biphasic composite with hierarchical structure. Right: 
viscoelastic behavior of the composite represented by a set of nested viscoelastic boxes. 


Working in the framework of the elementary theory of composite materials, we suppose that the 
connection between stress and strain can be obtained in one dimension, for this composite, using a 
sequence of nested viscoelastic boxes (Figure 1, right). The viscoelastic box of order n is composed of 
N,t+lelements in parallel: one linear viscous element with damping parameter 7, and N, complex 
viscoelastic elements. Each of these complex elements is the in series combination of a linear elastic 
element with compliance J, and a viscoelastic box of order n+1. Then, ifo,,(¢)is the stress applied to 


the viscoelastic box of order n at time?, and if ¢, (t) is its corresponding strain response, we have: 


é,(t)=[,0,(t-2)-o,,(¢’)-de'. (1) 


We assume that the stress and strain are always zero when t<0. The creep kernel of the viscoelastic 
boxes of the level n in the hierarchical viscoelastic structure is Q(t). The viscoelastic box of order 1, 


which contains all the other boxes nested in its interior, is supposed to be in series with a single elastic 
element of compliance /,. This element gives the impact or immediate elastic response of the 


considered element of volume of the composite material. If ¢", (t) is the strain produced in the elastic 


element of compliance J, , then: 


n? 


Enlt)=e'n()+ Enuld)- (2) 


But E', («)= Jn ‘One (0) » SO that: 


En (:) =In On (¢)+ [Que (¢ ~ t'): On (“’) -dt'. (3) 


But in parallel to the combination of elastic compliance with viscoelastic box of order n there is a 
viscous element. The stress acting through it verifies: 


Cw lt) =e = Fat). (4) 
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If¢, is the volume fraction of this viscous element, referred to the whole volume of box number n, 


and if g, is the volume fraction of the in series combination of an elastic element with box number 


n+1, from the elementary theory of composite materials we have, applying the popular rule of 
mixtures: 


o,(()=N,, “Pn “Ons (t)+ Pry “Ony(t). (5) 
The volume fractions verify, for every n: 
Ny ‘Ont ny =I. (6) 


The possibility of applying the rule of mixtures, at least as a convenient approximation, can be 
justified as follows. If A, is a cross sectional area corresponding to the boxes of level n in the 


viscoelastic hierarchy and A,, is the area corresponding to the dashpot, then: 


A, : on (r) = Ni, i Ans 7 Ons (1)+ Aas i o,,,(t) SO On (t) = N,, . An : Ont (1) | — . ony (t) * 


nv 


A 
so we have N,,-—#+ 
A, A, 


certain component in a set of parallel cross sections of a composite body can be used as an estimation 
of the volume fraction of the component in the body [9]. As consequence of this connection between 


AS A, =N,,° Ani + 4 


n+l nv 9 


=1 The average fraction of the area occupied by a 


area averages and volumes of a given phase in a body of composite material, if he is interpreted as 


n 


an estimation of the average fraction of parallel cross sectional areas of viscoelastic boxes of level 


n+1 inside a box of level, we have 9, = aa An analogous derivation gives 9,,, = E Of course 
both stresses and strains must be interpreted also as suitable averages. 
A Transmission Line Equation for the Viscoelastic Hierarchy 
From (5) we have: 
Pny -Ony(t)=o,(t)—N, “Py “On (t) © (7) 


From (23), taking into account equations (2) and (4), «,,(t)=«"', (()+¢,,(0) ando,, (=n, Ze, (t), as 


well as the constitutive relation between stress and strain for the in series spring inside boxes of level 
n, that ise’, (t)=/,,-,,,(t) , we obtain the following system of linear differential equations: 


0 ()= (<,()—«,(¢)) 


En 
Ot J 


Gol-erald) (8) 


Diy Ws Na Op 


n-l n 


The idea now is to approximate (8) by a partial differential equation of the transmission line type 
used in electric engineering, following a general method and employing the calculations developed in 
reference [10]. 

We take the Laplace transform in both members of (8): 

(s) = G4 (s)- En (s)) é,, (s)- Ena (s)) 2 (9) 


—N,, “Dy * 
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After some reordering, (9) may be put as follows, being 7, =@9,,-7,-J, : 


J 


n 


Ty 8-B(0)=Byal0)-2-Fy() Eyal o)o{ 1-2 Be Ze) (66) —256) 


If the branching numbers, volume fractions and elastic compliances are constants then: 


T,,-8-E,(8)=&,4(s)—2-&,(s)+ E,4(s)+ -N-)- E,(s)-E,uls)) - (10) 
Following the approach of [10], let us define ¢ = —— and introduce a new variable: 
s-() =e (11) 
lye) te” 


If the parameters of the dashpots grow like 77, [5] -7, taking into account that, from (11), 
P 


n= ine we obtain: 
Ing 


1 
Ke nln} — nx 
5 =e (a) =e Ing =x, (12) 
Here, by definition: 


gr F 
v= tn(v-o'*) ‘ (13) 
1 
——— 
{y a 
As we shall see in the discussion, v is always positive. If we take nas a continuous real variable, 
then x= &" is a regular real function ofn. We define: @,(s)=€(x,s)= a(e",s) 


Then 5(0)-B4(0)=8{ 1] ass) may be approximated by -=-(é 2a Es) : 


é,(s)—&,,,(s) = @(x,s)-&(E-x,s) may be approximated by — x-(€-1)- < &(x,s)and @, ,(s)—2-é,(s)+é,,,(s) 
x 


n+l 
may be approximated by x-(é-1) E(x 23 («s)] . Taking all this into account, after several 
IX IX 
operations we obtain the following transmission line equation in Laplace’s domain: 
Os .e 
&(x,s)—7, -s-x" 1 -E(x,5)x0. (14) 


ax? 


Here, with é = pied 
N-@ 


; . 15 
* (€-1)7 -€74 = 
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Equation (14) allows an approximate description of a finite viscoelastic hierarchy, say of N 
levels. This differential equation is equivalent to a modified Bessel equation and from a mathematical 
point of view is identical with the self-similar transmission line equation derived in [10]. The 
solutions of equation (14) are linear combinations of modified Bessel functions, both of the first kind 


a 1 aoe 
x fants” and of the second kind x” Afni } The parameter _p verifies: 


1 
fe 
jot, (16) 
l+v 1 
Hf 

Q 


When rf, -s is positive and much less than 1, the dominant term in the formula of the Laplace’s 


Power Law Strain 


transform Q,(s) for the viscoelastic creep kernel Q,(r) can be estimated by a suitable modification of 


the corresponding formula for the input impedance of an infinite self-similar transmission line that 
appears as case C in the mathematical appendix to chapter 7 in reference [10]. We obtain: 


O(s)es-prrr. Te) ot (17) 


Here, forx >0, [(x)=fe™“-w*'-du is the Gamma Function. Now we don’t have indeterminate 
0 


parameters as we had as a result of applying the method of functional equations [8]. 


Stretched Exponential Stress Relaxation 


If a linear viscoelastic material has both instantaneous (elastic) ¢,(t) and delayed ¢,(¢) strain 
response to an applied stress o(t) that is zero for negative times, the total strain ¢, (t)+ Ey (t) is given 
by: 

e(t)=Jo-o(t)+ [K(t-1')-o(t')-dt'. (18) 


If the case of the composite material studied in this work, the creep kernel K(t) may be identified 
with Q,(t), the kernel of the box in level one in the viscoelastic hierarchy. Inverting (18) and 
introducing an elastic impact modulus M, such that/,-M, =1 and a stress relaxation kernel L(t): 


o(t)= My -e(t)-[.L(t-2*)-e(¢’)- ar’. (19) 


Taking Laplace's transforms in (18) and (19), it is possible to relate the transform Z(s) of the stress 
relaxation kernel with the transform K(s) of the creep kernel: 


(20) 
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Now we substitute the Laplace's transform of the power law kernel (see equation (17)) in place of 
K(s) in (20) and we obtain: 


Es) | al , (21) 


12) 0 r(p) lg 22) 


Inverting the transform (21): 


i es { i)" al zt) (23) 


As shown in [11], the function Sw ” can be approximated by p-(2- p)? -e? 2-7" 


ol ((n + 1) Pp 


Then: 


(24) 


l-p 
= My Z{2) -p.-p) 
Vans 

From (24) the relaxation of stress at a constant strain ¢, for a viscoelastic body is given by 
o(t)=M(t)-e) with, by defintion: M(r)= (v1, -(‘1(+)-at', Taking into account (24) it follows that the 


stress relaxes according to a stretched exponential: 


o(t)=M(t)-€) =o a. : (25) 


In this equation: o) =M,)-€) A formula for the characteristic stress relaxation time 7, given as 
function of the characteristic time t, of the hyperbolic strain response follows from formulae (22) 
and (24): 

= fs , (26) 


The parameter p of the power law strain response. 


From (16) we conclude that p must be positive, because both & = <a and + are greater than 1. 
. QP — 


But from both, experimental results and thermodynamic requirements related with elastic energy 
dissipation, we know that p must be less than 1. This is equivalent to the restriction v>1. Taking 


Advanced Materials Research Vol. 853 15 


into account the definition of v given in equation (13), as well as the inequality V-@<1, we see 
that v>1 is verified if and only ifN-g'%>1. The physical origin of this last inequality and the 
physical origin of the scale law for the dashpot viscous modulus 77, (9° ie ~n,, can be found 


examining the geometric scaling laws of the viscoelastic hierarchy. Let us consider an element of 
volume V, in the composite material corresponding to a combination in series of an elastic element 


and a whole viscoelastic box of level n. Let us introduce a characteristic length/, corresponding to 


volume V, and a characteristic area S$, ~/? of the interfaces between the volume element and its 


. A 2 | 2 
viscous matrix environment. We assume that /,,, ~ oe 1,5 SO Siyy * of -S, and V,,,~@-V, . 
Furthermore, we suppose that the force F, taken by the viscous matrix in the volume element is given 


O€ 
by F, =u-S,-— 
Yr, fad n ot 


. Then, the force per unit volume a, is given by: 


Fi a ay". Ss Oe 
Gy =a =e = Sho I 3 He = . We suppose now that the viscous stress o7,,is related 
n t 1 t 


0E,, 


with the density of volumetric viscous forces by a constant factor: o,,=y-6,, Sothato,,=7, ave 


1 


n-1 
with 77, {| -‘m and m=y: wo . As consequence, we obtain the scale law of viscous 
1 


dissipation 7, (o8 7 =n with B= : . According to this model, the parameter of the power law and 


ire] 


of the stretched exponential verifies p=—+——<. Now v is positive and p is less than one if 


file? 


; 2 ’ ; pens 
N Sat ce N go =(v gh ) is greater than 1: this happens if the total area related with viscous 


dissipation increases from one level of the hierarchy to the next one. 


Discussion and Conclusions 


A unified phenomenological approach to power law strain and stretched exponential stress relaxation 
was developed in this paper and in reference [8], for a simplified example of viscoelastic hierarchy. 
An analytical formula for the power law parameter p as function of the branching number and 


volume fractions of the hierarchical composite material was derived for a simplified example. In this 
case the parameter of the stretched exponential turned to be the same as the power law parameter. 
The characteristic times 7, and t, were given as well defined functions of the mechanical 


parameters of the composite material. Due to its scale invariance, the hyperbolic strain response has 

an infinite and unbounded set of local time scales _. The characteristic time T, must not be 
P 

considered as a time scale of the hyperbolic response. It gives a measure of the numerical order of the 


interval of time (since the application of a step in stress) after which the power law can be applied. 
The characteristic time 7, gives a measure of the numerical order of the interval of time (since the 


application of a step in strain) after which the stretched exponential law for stress relaxation can be 
applied. The stretched exponential also has an infinite and unbounded set of local time scales 


Pp 
(*f } -t Contrary to the hyperbolic response case, now t, appears in formula for the local time 


scale of stress relaxation [8]. 
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If properly validated in practice, the kind of approach developed here could be used as a guide to 
the design of power law viscoelastic materials with desired values of the power law parameter p and 
characteristic timesz, andr, . 

The transmission line approach enables us to model different in series and in parallel viscoelastic 
hierarchies, considered as generalized quadrupoles. In this case strain corresponds to voltage and 
mechanical force (stress multiplied by area) corresponds to current. Dashpots correspond to 
capacitors and elastic elements correspond to resistances. With these correspondences, the 
well-developed mathematical tools of distributed parameter electric networks can be applied to study 
viscoelastic hierarchies. Random branching number and moduli could be introduced. 

Finite viscoelastic hierarchies show, besides the high frequency cut off that corresponds to the 
beginning of the deformation process, a low frequency cut off for long times of deformation. In the 
model based in a set of nested boxes, the low frequency cut off behaviour depends of the termination 
(elastic elements or dashpots) of the last box in the hierarchy. 
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